Asymptotic boundary layer method for unstable trajectories: Semiclassics for 

individual scar wavefunctions 
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We extend the asymptotic boundary layer (ABL) method, originally developed for stable resonator 
modes, to the description of individual wavefunctions localized around unstable periodic orbits. The 
formalism applies to the description of scar states in fully or partially chaotic quantum systems, and 
also allows for the presence of smooth and sharp potentials, as well as magnetic fields. We argue that 
the separatrix wave function provides the largest contribution to the scars on a single wave function. 
This agrees with earlier results on the wave-function asymptotics and on the quantization condition 
of the scar states. Predictions of the ABL formalism are compared with the exact numerical solution 
for a strip resonator with a parabolic confinement potential and a magnetic field. 

PACS numbers: 03.65.Sq, 05.45.Mt, 73.21.-b 



INTRODUCTION 



Semiclassical methods provide valuable insight into 
the eigenmodes of electromagnetic microresonators and 
mesoscopic electronic devices because they establish a 
direct relation to classical trajectories. For stable or- 
bits the corresponding wavefunctions are localized in the 
vicinity of the trajectory and can be constructed by the 
elegant asymptotic boundary layer (ABL) method, which 
is based on rectified Gaussian beams and their harmonic 
transverse excitations p], 0, 0, 0] . Describing wavefunc- 
tions associated with unstable orbits is a considerably 
more difficult problem. For classically chaotic quantum 
systems [5[ , the rapid stretching and folding of the phase 
space translates into a randomization of most eigenstates, 
which can be addressed through statistical descriptions 
0]. However, some eigenstates deviate significantly 
from a random wave pattern in that they are enhanced 
in the vicinity of unstable orbits [l, 0, 0] . These so-called 
scars are most pronounced for short periodic orbits, and 
often dominate the individual properties of quantum sys- 
tems. Their influence has been detected in numerous nu- 
merical investigations and experiments on a large variety 
of systems, such as microwave cavities [10[, resonant tun- 
nelling diodes [11], and, in particular, in the lasing modes 
of dielectric microresobators [121 ]. 

Pioneering theoretical investigations of individual scar 
wavefunctions retained the assumption of a Gaussian 
transverse profile and its harmonic excitations @, 0, Ojjl, 
[HI . It was soon realized that this requires to superpose a 
large number of transverse excitations [HI, [l6j] . On aver- 
age, the large-distance asymptotic decay of the transverse 
profile is expected to follow a power law 
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Studies of hyperbolic fixed points (with effective period 
0) also point into the direction that the transverse pro- 
files are essentially non-Gaussian [2] 1 [2 2I (23[. On the 
other hand, numerical investigations [24fnave led to the 



conjecture that scar quantization only involves a single 
(longitudinal) quantum number, which means that the 
transverse mode profile is fixed. 

In this paper we clarify the nature of scarred states 
from the perspective the ABL method, which offers a sys- 
tematic semiclassical expansion for the Schrodinger equa- 
tion in the asymptotic boundary layer around a classical 
trajectory |l|. Compared to other semiclassical meth- 
ods, a main advantage of the ABL method is its appli- 
cability to a broad class systems, including systems in 
external potentials and magnetic fields. ABL expansion 
up to quadratic terms of the effective potential yields 
the Gaussian wave packet solution, originally proposed 
by Heller Q as a convenient ansatz to study dynamics 
of complex molecules. Subsequent mathematical studies 
P, 0, [30[ advanced the method into a powerful general 
tool to obtain asymptotic solutions to various equations 
of theoretical and mathematical physics, ranging from 
optics to gravity. 

Most of these applications, however, focus on station- 
ary eigenstates associated with stable periodic orbits. 
The application of the ABL formalism to unstable orbits 
faces the problem that Gaussian beams fail to satisfy the 
periodicity condition along such orbits. In principle, this 
problem can be circumvented by numerical superposition 
of a large (ideally, infinite) number of Gaussian solutions 
for an artificially stabilized orbit [15]. However, as we 
will show, it is possible to take a more direct analytic ap- 
proach, which focusses on the mathematical difference of 
the ABL equations for stable and unstable trajectories. 
This allows to obtain closed analytical expressions for in- 
dividual wavefunctions localized around unstable trajec- 
tories. A detailed analysis of these solutions reveals the 
special role of the so-called separatrix solution in the for- 
mation of scar states. This yields a transverse mode pro- 
file that corresponds to an appropriately rectified specific 
solution of the hyperbolic fixed-point problem [2l], [22|, [23[ 
and also satisfies the correct large-distance asymptotics 
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[I?], EH • The periodicity condition for the separatrix so- 
lution leads to a Bohr-Sommerfeld quantization formulae 
with a single longitudinal quantum number, in agreement 
with the earlier numerical conjecture of Ref. [24[. These 
general considerations are presented in Sections HHHV1 

In Section [V] we illustrate the flexibility of the ABL 
approach by applying it to the example of a two- 
dimensional strip resonator in a magnetic field and a 
parabolic confinement potential. Comparison of the ABL 
predictions with the results of numerical computations 
demonstrates excellent accuracy of the method. 

Section lVll contains discussion and conclusions. Details 
of the derivations are presented in the Appendix. 



II. ABL FORMALISM 

In the following two sections we provide the main re- 
sults of the ABL formalism and cast them into an outfit 
which naturally contrasts states localized around stable 
and unstable trajectories. Details of the derivation of 
the main formulas are given in Appendix. To make the 
presentation compact we concentrate on two-dimensional 
quantum systems in a constant perpendicular magnetic 
field B, with vector potential A = B{—y,x)/2. 

We employ a curvilinear system of coordinates asso- 
ciated with the main classical trajectory tq(s), which is 
parameterized by its arclength s. A pair of coordinates 
(s,n) defines the position vector in the vicinity of the 
trajectory as r = r (s) + e n (s)n, where e n is the normal 
unit vector and n measures the distance from the trajec- 
tory. The ABL formalism seeks a solution in an asymp- 
totically small boundary layer of width oc around the 
main trajectory, which introduces a natural scaling of the 
distance as n = \ffiv. The solution of the Schrodinger 
equation is sought in the form 

00 

v) = e ^ s oW+Si(s)VKu) J- h j ^ j \s : v\ (1) 
3=0 

which represents an expansion in half-integer powers of 
Planck's constant h (in contrast to the text-book WKB 
expansion in integer powers of Ti). This expansion of the 
Schrodinger equation (see Appendix) yields the exponent 
in Eq. (pQ) in the form of the classical action functions 

50 = f {a(s) + e(A x el - A y e x n )} ds, (2a) 

5 1 = e(A x e x n + A ye y n ), (2b) 
a(s) = y/2m[E-u(r (s))], (2c) 

where u(r) is the scalar potential. For one obtains 
a hierarchical set of equations. If the expansion of the 
effective potential in the system can be restricted to the 
terms of second order in n (the standard assumption in 
the linear theories [9|, LL4[ ) , one can restrict the series in 
Eq. (pQ) to the leading term ip = and obtains the 



so-called boundary layer Schrodinger equation 

^(-,")=(~« + Y^-|a)^(-."). (3) 

7/ N 2mu2 m 2 u 2 2mu\ eB 

d(s) = — — + — . 

a z a 4 pa z pa 

Here the dot denotes the partial derivative <9 S , while p is 
the geometric radius of curvature of the trajectory. The 
functions are obtained from the expansion 

u(r) w u (r (s)) + ui(s)Vhv + U2(s)Tiv 2 . (4) 

Equation (j3j) is analogous to the Schrodinger equation 
for a one-dimensional nonstationary oscillator, where s, 
a(s), and d(s) take the role of time, mass, and harmonic 
frequency, respectively. As shown in the Appendix, a 
general form of its partial solutions can be constructed by 
establishing their relation with the classical trajectories 
in the vicinity of the main one. The final result yields a 
set of solutions of the form 

«->«^^(^)^ 

(5) 

where are parabolic cylinder functions with arbitrary 
index £, while T = p/z and T = p/z are defined by 
two independent solutions of the Hamilton equations in 
variation 

z = -, p— —adz, (6) 
a 

that satisfy pz — pz = w (w is the Wronskian, which here 
is a constant). A pair (z,p) describes the classical trajec- 
tories in the vicinity of the main trajectory. Once Eqs. 
(|6j) are solved, Eqs. ([1]), ([2]) and ([5]) define semiclassical 
solutions of the Schrodinger equation. A second class of 
partial solutions is found by interchanging z <-» z in Eq. 

The similarity of Eq. ([5j) with a standard oscillator 
can be seen by assuming constant a ("mass" m) and d 
(squared "frequency" uo 2 ). If one also chooses solutions to 
Eq. (|6]) as z oc exp(zcjs), Eq. ([5]) yields familiar oscillator 
wave functions. In the general case, the detailed depen- 
dence of a and d on s results in nontrivial modulations 
of the wave function profile along the trajectory. 

III. PERIODIC ORBITS 

Since there are two partial solutions of Eq. (|6|), the 
general solution Eq. ([5]) for the wave function possesses 
one degree of freedom. For closed periodic orbits, how- 
ever, this freedom is restricted by an additional period- 
icity condition. As we shall see, this condition naturally 
distinguishes between stable and unstable orbits. 

For stable orbits we follow a standard procedure, the 
first step of which is to find the Floquet solutions of Eq. 
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(j6]), which fulfill the conditions z(T) = exp(z0)z(O) and 
p(T) = exp(i0)p(O). Here T is the period of the orbit 
and is a real Floquet phase. If (^1,2,^1,2) are two inde- 
pendent solutions with initial conditions (1,0) and (0, 1), 
respectively, the general solution is a linear combination 
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(7) 



Floquet solutions are found by determining ai, a2, and 
A = exp(i(/>) from the eigenvalue problem 



M 
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Zl (T) z 2 {T) 
Pi(T) P2 (T) 



■ (8) 



Here M is the monodromy matrix (with a unit determi- 
nant), which depends on the stability of the orbit. 

For stable orbits |Tr M\ < 2, and M has two com- 
plex conjugate eigenvalues with a unit absolute value. 
The corresponding Floquet solutions are also complex 
conjugate, (z,p) and (z,p) = (z*,p*). We use the con- 
vention that (z,p) is the solution with ImT > 0. This 
choice is well-defined since ImT ex \z\~ 2 does not change 
sign along the trajectory. With this choice, Eq. ([5]) is 
restricted to normalizable harmonic oscillator functions, 
which are also referred to as the Gaussian beam (or ray) 
solutions p], [ij] . The periodicity condition for the solution 
defined by Eqs. (pQ), ([2]) and ([5]) is obtained from the phase 
increment for a single round trip along the orbit, which 
depends on the classical action So and the phases of the 
various powers of z,z in the expressions. This yields the 
Bohr-Sommerfeld quantization formula for stable orbits, 
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a(s)ds + $ = Ti[2im + (m + 1/2) 0(E)], (9) 



where <£> is the magnetic flux through the orbit, while n 
and m are longitudinal and transverse quantum numbers. 

For unstable orbits, |Tr M\ > 2. In this case the 
eigenvectors and eigenvalues of M are real, such that 
A± = exp(±A) where A is the dimensionless Lyapunov 
exponent (corresponding to a time-domain Lyapunov ex- 
ponent At = A/T). For this case the harmonic oscillator 
functions with integer index m in Eq. ([5]) violate the pe- 
riodicity condition. However, the periodicity can still be 
achieved with the choice of the £ = — 1/2 + ir] (where 77 
is real). The resulting functions are similar to those for 
an inverted (negative) harmonic potential. 

The periodicity requirement of the wave function then 
results in the Bohr-Sommerfeld like quantization condi- 
tion for unstable trajectories 



/ a(s)ds + $ = ft [27m ± rj\(E) + na/2]. 
Jo 



(10) 



The Maslov index a counts singular points along the tra- 
jectory, where z or z vanishes. Since T oc these 
points produce a square-integrable singularity oc z -1 / 4 in 
the wave function [see also Eq. ([T3]h below], which can 
be further regularized by uniform approximations (see, 
e.g., Refs. lallsl)- 



IV. CONTRIBUTION TO SCAR STATES 

The ABL solutions defined by Eqs. ([5]) and (|6]) are only 
valid in a boundary layer of width \n\ = y/h\u\ = 0(Vh) 
around the orbit. In order to construct eigenstates as- 
sociated with the unstable trajectory, the semiclassical 
wave function has to be matched to a quasi-random back- 
ground beyond this layer [26[ , which in principle can mix 
solutions with different values of 77 and n. The number 
of contributing solutions is defined by the interplay be- 
tween the coupling of the solutions with the background, 
and the energetics of the different states as determined 
by the quantization condition in Eq. ([TO]) . In essence, one 
therefore deals with a scattering problem. As confirmed 
in recent works on quantum resonance wave functions 
OS H3 5 HH, an important scale in this context is the 
Ehrenfest time [29[ , defined by the time of a wavepacket 
of initial size of de Broglie wave length to spread across 
the entire accessible phase space of the system, 



tEhr = A T In K 



P h. 



(11) 



Here At = A/T is the time-domain Lyapunov exponent, 
and TVph = L\j2mE/l% ^> 1 is a dimensionless measure 
(in units of Planck's constant) of the total accessible vol- 
ume of the phase space, where L is the characteristic size 
of the system. 

The Ehrenfest time defines an energy window 



E = E no ± 



27rh 



(12) 



over which partial ABL solutions are strongly mixed 
among each other via the background states. Here 
E no — E nQy7] =o the characteristic energy following from 
the Bohr-Sommerfeld quantization condition (jTUJ) , with 
r] = 0. Earlier numerical observations provide evidence 
that scars are quantized with energies close to E no (24[. 
This is the energy of the symmetric separatrix solution, 



-<-''^V^- 7 1 ' U?)'" ' ' (13) 



Separatrix solutions without trajectory-specific z and z 
were earlier discussed in the context of scarring phe- 
nomenon as the solution of the hyperbolic fixed-point 
problem [2l|, S El ■ The asymptotic form of Eq. (O 
for large v 



1 fwv 2 \ 



(14) 



recovers the well known result of Bogomolny [17[ (again, 
generalized to contain the trajectory-specific z\ which is 
often interpreted as the profile of a scarred wave func- 
tion associated with few closely located unstable orbits 
We also note that the description of scars in earlier 
works [15[ employed the ABL formalism for stable orbits 
together with the variational procedure to minimize the 
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squared transverse energy. This effectively determines 
the wave function closest to the separatrix solution ([T3]h 
but represents this as a sum of harmonic oscillator func- 
tions. 

Expanding the quantization condition (fT0|) around 



E = E no one obtains 



EK,E no + —8n + hX T V, 



(15) 



where Sn = n — uq. The approximate 'level spacing' of 
states with different rj can be estimated from the asymp- 
totics of Eq. ([5]) at v — > oo, 



/ / \ i I r 2 . / /r-r 

^ ^) oc — p= exp < z — ^ —277m I // \ 



(16) 



Assuming that the mixing occurs when the state spreads 
to a transverse distance L tr = O(L), the typical spacing 
of solutions with different 77 is 
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(17) 



probability from the orbit boundary layer. For unstable 
orbits the localization is given by a power law, follow- 
ing from the asymptotic expansion in Eq. ([16|) ; the ABL 
wave functions then are not square integrable, and the 
escape rate increases logarithmically with the system di- 
mensions. The symmetric separatrix solution ([T3|) com- 
bines a large maximal value in the vicinity of the orbit 
with a small transverse probability current. At 77 < 
the maximum of Eq. ([5]) shifts from the orbit, eventually 
moving out of the boundary layer. At 77 > the trans- 
verse profile of the solution becomes flat and its weight 
shifts from the main trajectory. [Interestingly, the sepa- 
ratrix solution is also best behaved in the vicinity of focal 
points (where zz = 0); for 77 7^ the ABL solutions dis- 
play a jump by a factor exp(?77r).] The dominance of the 
separatrix solution therefore has a simple physical ori- 
gin: the 77-dependent term in the quantization condition 
([TO]) controls the transverse momentum of the solutions, 
which in turn determines the leakage out of the boundary 
layer. 



V. APPLICATION TO A QUANTUM 
RESONATOR MODEL 



where () s denotes averaging over the trajectory. 

Under certain conditions, only a single solution ([5]) is 
found in the Ehrenfest window ([T2]) . From the longitu- 
dinal quantization one obtains the standard condition 



T < ^Ehr- 

The transverse quantization gap ([TT]) implies 



T< 0(X T 



(18) 



(19) 



with C > depending on the precise choice of L tr . This 
condition is satisfied for all T < tehr, and therefore 
is weaker than the longitudinal quantization condition. 
Consequently, as long as the period of the orbit T is less 
than the Ehrenfest time, only a single ABL solution will 
contribute to the semiclassical eigenstate. 

The arguments from the above procedure of matching 
of the ABL solution to the random background therefore 
recovers the main phenomenological result on the strong 
scarring along unstable trajectories, Eq. (fTSjl . which was 
earlier obtained by analysis of time-dependent Gaussian 
packets [34[ . Because of the conceptual relation to a scat- 
tering problem, this argumentation can be enforced by 
taking the transverse probability flux across the scar into 
account. This can be studied using the probability cur- 
rent in curvilinear coordinates, 



J = UsJu) = ( ail>*ip, - {^ v d^* - ^*d^ v } ) , (20a) 



d s j s + dv3v = 0. 



(20b) 



For truly stable orbits the wave functions are expo- 
nentially localized, resulting in a minimal leakage of the 



In order to assess the predictive power of the above 
ABL formalism for the description of individual scarred 
eigenfunctions we now turn to a specific system: a two- 
dimensional strip resonator (quantum wire) in a longi- 
tudinal confinement potential and a perpendicular mag- 
netic field B. In the x direction, the system is confined by 
two parallel impenetrable walls at x = and x = d, while 
in the y direction the confining potential is parabolic, 
U(y) = mcjQy 2 /2. In the computations we set d = 10 Ib 
and c^o = 2c^ c , where uj c = \eB\/m and Ib = \/Ti/muu c 
are the cyclotron frequency and the Landau magnetic 
length, respectively. 

The classical dynamics of this system is of the generic, 
mixed, type, where stable orbits coexist with unstable 
orbits [see phase space portraits in the middle panels of 
Figs. Ufa)- (d)]. It was demonstrated previously for a sim- 
ilar system [4J that the ABL method provides an accurate 
description of quantum states associated with stable or- 
bits. Here we analyze an energy-dependent family of un- 
stable bell-shaped orbits with two reflection points (0, 0) 
and (d, 0) at the hard walls [see white curves in the right 
panels of Figs. [2(a)- (d)]. The orbits consist of an upper 
and a lower arc, which are symmetrical to each other and 
are indexed in the following by i = 1,2. 

We solve the equation in variation ([6]) on each arc and 
then link these solutions together using a reflection ma- 
trix, which for non-vanishing magnetic field takes the 
form 



R 



-1 







-2a; c tan0 — 1 



(21) 



where is the angle of reflection. The monodromy ma- 
trix of the full orbit is found from the product M = 



5 



( a ) E=90.754fico c E=90.795fio) c ( b ) E=92.474ficD c E=92.554R« c 
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FIG. 1: (Color), (a)-(d): Sequence of scarred wavefunctions localized around the bell-shaped unstable orbit in an electronic 
strip resonator with a parabolic longitudinal potential and a transverse magnetic field. Left panels: Color-coded modulus of 
the exact wavefunctions \^(x,y)\ (blue: small amplitude, red: large amplitude). Right panels: The corresponding semiclassical 
ABL wavefunctions. The white curves depict the bell-shaped orbit, computed at the semiclassical energy. Upper middle panels: 
color-coded Husimi representations of the exact wavefunctions. Lower middle panels: classical phase space portraits, computed 
at the energy of the exact wavefunctions. The red dot denotes the position of the bell-shaped orbit. The white circle in the 
Husimi representations delimits the energetically accessible region of the classical phase space. 



RM.2RM\^ where Mi is the fundamental matrix of arc 
i. There are eight singular points, where zz = 0: four 
well isolated focal points around x « (0.5 ± 0.15) 
y « ±0.5 <i, and two pairs of closely spaced focal points 
around x « (0.5 ± 0.01) d, y « ±0.65 d. The stretching 
factor for the orbit (the largest eigenvalue of the mon- 
odromy matrix) is A w 4.6, corresponding to A « 1.53. 

Semiclassically, the eigenstate is found as the sum of 
the two separatrix solutions ([T3]) along each arc, 

* = s, v) ± r^ 2 (^2, s, v). (22) 

Here the coefficient r follows from the continuity condi- 
tion at the reflection points, and simply takes the value 
r = — 1 for unstable orbits. The condition \£ = at 
the boundaries x = and x = d of the resonator can 
be taken into account by adopting a mirror reflection 
method [3, 0]. In this method, the wave function with 
the required boundary condition is obtained as the dif- 
ference between the original solution and its mirror re- 
flection at the boundary. In principle, the presence of the 
magnetic field requires to introduce an additional gauge 
field to eliminate the phase jump associated to the differ- 



ence T2 — Ti = 2cu c entering the off-diagonal element of 
the reflection matrix Eq. ([2T]) . However, this affects only 
a small neighborhood ex Vh of the point of reflection. 

In order to compare the resulting ABL solutions with 
exact quantum mechanics, we construct numerical eigen- 
states in an orthogonalized basis spanned by the exact so- 
lutions of the quasi-one dimensional system without the 
hard walls, which separates in the gauge A = —Bye x . 
The boundary conditions ^(x = 0, d; y) — are imposed 
via a singular value decomposition (for details on such 
methods see, e.g., Ref. |33|). 

Our numerical computations reveal that the family of 
bell-shaped unstable orbits supports a long sequence of 
scarred wavefunctions, which are found at almost equidis- 
tantly spaced energies. Four consecutive examples (with 
longitudinal quantum numbers n = 66, 67, 68, 69) are 
shown in the left panels of Figs. [2(a)- (d). The corre- 
sponding ABL wavefunctions (right panels) accurately 
capture the typical spatial extent of the scar signature, 
including the position of the focal points, where the scar 
structure shrinks while the amplitude is significantly en- 
hanced (close to the almost-degenerate pairs of focal 
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(d) 




FIG. 2: Profile of the exact and semiclassical wave functions 
calculated along the line y = 0. The width of the boundary 
layer where the ABL approach is valid is indicated in the 
upper left corner. 



points the wave function could be further regularized by 
adapting the general techniques of uniform approxima- 
tions; see, e.g., Refs. [3l|, |32|). 

The clear correspondence between the exact quan- 
tum states and the ABL separatrix solutions is fur- 
ther confirmed by the Husimi representations H(y,p y ) — 
\(y,p y \d x ty(x = 0,7/))| 2 , which are obtained by over- 
lapping the derivative of the exact wavefunction at the 
left wall with coherent states (minimal-uncertainty wave 
packets) that are parameterized in Birkhoff coordinates 
y,p y . The Husimi representation, shown in the upper 
middle panels of Figs. [T](a)-(d), in all cases displays a 
clear maximum at the position of the bell-shaped orbit 
(red dot in the chaotic part of the phase-space portraits), 
which lies in the chaotic part of phase space. 

Besides the visual agreement, the separatrix solutions 
also recover the following essential characteristics of the 
numerical solutions: (i) The wavefunctions are of the 
right symmetry with respect to inversion around x = 
0.5 y = (as there are eight focal points, wavefunc- 
tions are symmetric under this operation when n is odd, 
while they are antisymmetric when n is even), (ii) The 
profile of the wave function along the line y = extracted 
from the ABL and the numerical solution (after spatial 
averaging over a range of 0.05 d to eliminate the speckle 
fluctuations due to the random background in the numer- 
ical solution) shows a very good agreement within valid- 
ity region of the ABL solution, x/d < 0.15. The ABL 
captures well both the period and the relative strength 
of the oscillations in the numerical solution, (iii) The 
exact energies are in excellent agreement with E no , i.e., 
the semiclassical prediction from Eq. ([TO)) with r] = 0. 
These energies are almost equidistantly spaced, and the 
semiclassical error is less than 5% of this spacing (the 
error is also small compared to the mean level spacing 
of all states), (iv) For each longitudinal wavenumber n 
the numerical computations only deliver a single scarred 



wavefunction localized around the unstable trajectory (in 
contrast, many transversely excited wavefunctions are 
supported by the stable trajectories in the island around 
y « 0, p y « 0.3muod). In other words, we do not find 
any other scarred states on this orbit. Therefore, the 
actual scar quantization window is reduced beyond the 
constraint obtained by the propagation of time depen- 
dent Gaussian states [34], which predicts scarring in a 
larger window AE w hXr- However, one cannot exclude 
that signatures of off-resonant scarring show up in the 
wavefunction statistics of the states in this larger win- 
dow, which is a question beyond our focus on individual 
wave functions. 



VI. CONCLUSIONS 

In this work we extended the asymptotic bound- 
ary layer (ABL) method for semiclassical wave func- 
tions, originally developed for the description of Gaus- 
sian beams guided by stable trajectories, to the case of 
scar-like states localized in the vicinity of unstable peri- 
odic orbits. We focus on wavefunctions scarred by a sin- 
gle short trajectory and derive expressions of a universal 
form, valid up to order ft 3 / 2 , which apply to general sys- 
tems which may combine hard walls, external potentials 
and magnetic fields. The ABL equations are formulated 
in a curvilinear coordinate system associated with the 
classical trajectory. The system- and orbit-specific infor- 
mation enters via the solution of the classical equations 
in variations, which describe the stability of the trajec- 
tory. At fixed energy, the profile of the wavefunction is 
determined by a periodicity condition, which results in a 
Bohr-Sommerfeld quantization formula. 

Far away from the guiding trajectory, the ABL wave 
function decays as the inverse square-root of the distance, 
and therefore is not normalizable, but couples to the 
quasi-random background of other modes in the system. 
This coupling is the weakest for the symmetric separatrix 
state, which is associated to a particulary simple quan- 
tization condition. Since the separatrix state is also the 
most visibly localized ABL wavefunction, one can expect 
that it typically provides the dominant contribution to 
scars guided by a single short trajectory. In general, we 
estimate that the number of wave functions of different 
profile participating in the scar formation increases log- 
arithmically with the system dimensions, and therefore 
depends on the relation of the Ehrenfest time and the 
period of the trajectory. 

We verified our conclusions by comparison with exact 
numerical solution for a quantum strip resonator with the 
quadratic confining potential and a perpendicular mag- 
netic field. In this system, a family of unstable periodic 
orbits supports a long sequence of scarred states, with 
energies in excellent agreement with the simple quanti- 
zation condition for the separatrix solution. We find that 
this solution captures all the essential characteristics of 
the exact scarred states. The studied system is of the 
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generic type, with a mixed phase space, which raises the 
expectation that the ABL formalism is applicable to a 
large class of quantum- dynamical systems. 

An open question concerns the generalization of the 
formalism to scars supported by many trajectories, as 
well as the contribution of ABL solutions to off-resonant 
scars (typically revealed in wavefunction statistics). 

This work was supported by the European Commission 
via Marie-Curie excellence grant No. MEXT-CT-2005- 
023778. 



APPENDIX A: ABL EXPANSION FOR 
SCHRODINGER EQUATION 



for 5o(s), Si(s), and ^^(s, v). If the expansion of the ef- 
fective potential can limited to second order (a standard 
assumption in linear scar theories Q3] ) , one can termi- 
nate the series in Eq. ([T|) at the leading term ^^(s^v). 
This is equivalent to obtaining a semiclassical wave func- 
tion to the order of h 3 ^ 2 The first three terms of the 
series expansion of Eq. (|A3p read 



h~ 2 : A ^ (0) = 0, (A4a) 

h~ 3/2 : A + 2iB ^ 0) + Cv ^ (0) = 0, (A4b) 

h' 1 : Cv + 2iB ipW + V = 0, (A4c) 

with coefficients 



In this appendix we outline the main steps in the 
derivation of the ABL formalism. Further details on the 
formalism can be found in Refs. [3,[2[. Here, we present a 
straightforward derivation based on the direct expansion 
of the Schrodinger equation. 

The main trajectory ro(s), which is parameterized by 
its arclength s, generates an orthogonal coordinate sys- 
tem spanned by the normal unit vector e n (s) and the 
longitudinal unit vector e £ (s), which we express in the 
component form 



>(*) 



Recalling the conversion formulas for the partial deriva- 
tives in curvilinear coordinates, 

d x = y<9 s + j x d n , d y = -y<9 s + jydn, (A2) 



where p(s) is the curvature radius of the trajectory, we 
write the Schrodinger equation as 



^ ss np _ ^ n 



ieB 



-f ~ + ^ + ¥ [(r e n ) + n] * s 

ieB, . T e 2 B , , 0t 2m . 

+ — (r e,)*„ - ^j( r o + e nn) 2 * = ^(U - E)V, 

(A3) 

where U(s,n) is the potential inside the resonator. 

The central idea of the ABL approach is to search the 
solution to Eq. (|A3j) in an asymptotically small boundary 
layer of the width ex Vh in the vicinity of the main trajec- 
tory. A straightforward way to derive the corresponding 
semiclassical expansion for the wave function follows two 
steps: (i) scaling the normal variable as n = Vhv and 
(ii) expanding the resulting equation in powers of Vn, 
as written in Eq. ([T]), where each term tp^ is uniquely 
defined. 

We substitute Eq. ([T]) into Eq. (|A3[) . expand the re- 
sulting equation, and match the coefficients of the re- 
sulting series, which yields a hierarchical set of equations 



A = a 2 



r%e 2 B z 



eB(roe05i - S\ + eB(r e„)So - S 2 , 

(A5a) 

(A5b) 



B= ^?(roe,) + Si, 
C = -i ^ - 2# #i - eBSl - e55i(r e n ) 



eB . x ^ n a 2 
— (r ej)5i - 2mui 

9 P 



e 2 B 2 
4p 



e 2 B 2 



(Al) V = 2ia ^(°) + id ^ (0) + - da 2 v 2 ^ 



(o) 



■(r e n ), 
(A5c) 
(A5d) 



Equation (|A4a|) yields A = 0. Equation (|A4bj) is solved 
by setting B = 0, since C = is satisfied identically if 
A = B = 0, which yields the classical action solution 
^j. Finally Eq. (|A4cj) is equivalent to V = 0, which 
yields Eq. ((3]). After the substitution <j> = a" 1/2 ^ (0) this 
equation take the form of a nonstationary Schrodinger 
equation, which we write as 



Lcj) = < id, 



d 2 ad 9 1 ( 

— v > 6 = 

2a 2 j Y 



0. 



(A6) 



The fact that Eq. (|SJ) is a solution to Eq. ([3|) can be 
checked by a direct substitution. Here we present a direct 
derivation, which provides deeper insight into the nature 
of the ABL solution. First, we obtain a "ground state" 
solution to Eq. (|A6p by using the ansatz 



0o = p(s) exp 



(-) y , 



(A7) 



which is similar to the thawed Gaussians introduced by 
Heller [25[ to describe the time evolution of localized 
quantum wave packets. Substituting Eq. (|A7|) into Eq. 
(|A6|) yields equations for the unknown quantities T and 
/i, 



. r 2 

r + — + ad = 0, 

a 



/i = exp 



l 2i « 



(A8a) 
(A8b) 
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The Ricatti equation (|A8aft is solved by a standard sub- 
stitution T = p/z, after which one finds equations for 
p and z as given by Eq. (|6]). As we mentioned above, 
this Hamiltonian equations in variations define relative 
coordinates z and momentum p for classical trajectories 
in the vicinity of the guiding trajectory. Once Eq. (|6]) 
is solved, Eq. (|A8b|) delivers \i = z~ x l 2 . Equation (j6j) 
has two linearly independent solutions, (z,p) and (z,p), 
which define the constant Wronskian pz — pz = w. 

The classical origin of Eq. (|6|) can be established from 
the reduced classical action constructed for the trajec- 
tory n(s) in the vicinity of the guiding trajectory. In 
curvilinear coordinates this action reads 



Jo 



S = I x/2m{E - U(r)) { ( 1 - ^ | + fr > ds 



1/2 



eB 



(xy - yx) ds, 



(A9) 



where x = xq(s) + r y x (s)n(s) and y = yo(s) + r y y (s)n(s). 
Expanding Eq. (|A9P to second order in the small quantity 
n, one obtains 

S = S - [ e n ds+- [ (ah 2 - adn 2 ) ds, (A10) 
Jo ^ J 



where 



e(s) = eB 



a au\ 



P 2(E - u ) 



(All) 



The linear term in this expansion vanishes identically, 
e = 0, because n = is also a trajectory (the main tra- 
jectory). The linear coefficient in Eq. (|A10|) is related 
to Eq. (|A5c|) as C = 2as, which explains why C van- 
ishes. The Euler-Lagrange equation for the action (|A10|) 
is equivalent to the Hamiltonian system ([6]). 



A sequence of partial solutions to Eq. ([3]) can be con- 
structed with the help of creation and annihilation oper- 
ators p, E3 



(A12a) 
(A12b) 



A = \ — (-i z d v -p v) , 



A 1 " = J— (-i zd v -pv), 



which satisfy the following algebra 

[A,L] = [A\L] = 0, [A,At] = l. (A13) 

General solutions are then obtained from ip m = A^ m ipQ. 
This procedure leads to a recurrence relation, which is 
solved by wave functions of the explicit form 



(j>m(s,v) (X 



z m (r-r) m / 2 



Hrr 



T-r 



2% 



v e 2 



(A14) 

where H m are Her mite polynomials. 

The presented derivation allows for an unambiguous 
analytical continuation of the recurrence relation and its 
solutions to arbitrary indices m — > £. This yields Eq. 
([5]), where the parabolic cylinder functions take the 
place of the Hermite polynomials H m , and the coefficient 
a _1//2 is accounted for. Finally, we note that recalling the 
definition of T and T as well as that for the Wronskian, 
the solution can be compactly written as 



fa(s,v) (X 



e i(r+r> 




(A15) 
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